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A traveling wave parametric amplifier (TWPA) composed of a transmission line made up of a 
chain of coupled asymmetric superconducting quantum interference devices (SQUIDs) is proposed. 

The unique nature of this transmission line is that its nonlinearity can be tuned with an external 
magnetic flux and can even change sign. This feature of the transmission line can be used to perform 
phase matching in a degenerate four-wave mixing process which can be utilized for parametric 
amplification of a weak signal in the presence of a strong pump. Numerical simulations of the 
TWPA design have shown that with tuning, phase matching can be achieved and an exponential 
gain as a function of the transmission line length can be realized. The flexibility of the proposed 
design can realize: compact TWPAs with less than 211 unit cells, signal gains greater than 20 dB, 

3 dB bandwidth greater than 5.4 GHz, and saturation powers up to -98 dBm. This amplifier design 
is well suited for multiplexed readout of quantum circuits or astronomical detectors in a compact 
configuration which can foster on-chip implementations. 

PACS numbers: 85.25.Am, 85.25.Cp, 84.30.Le 


I. INTRODUCTION 

Over the past decade Josephson parametric amplifiers 
have proven essential in experiments studying quantum 
jumpsi , tracking quantum trajectories^Ti^ and real-time 
monitoring and feed-back control^i^i^ of quantum bits 
(qubits). In parametric amplifiers high gain is achieved 
when the signal to be amplified interacts with a non¬ 
linear medium for as long as possible. State-of-the-art 
Josephson parametric amplifiers utilize resonant circuits 
to increase the interaction time of the signal with the non¬ 
linear medium2r— . As a consequence the instantaneous 
bandwidth and maximum input power allowed are sig¬ 
nificantly reduced. Such limitations have renewed inter¬ 
est in superconducting traveling wave parametric ampli¬ 
fiers (TWPA)iiA2 which achieve long interaction times 
with the nonlinear medium by extending the electrical 
length over a long transmission line instead of multiple 
bounces in a resonant circuit, as a result TWPAs do not 
suffer from the same bandwidth and dynamic range lim¬ 
itations that cavity based amplifiers do. The main chal¬ 
lenges in TWPA designs is that optimum gain is achieved 
when phase matching conditions are met. Supercon¬ 
ducting TWPAs have been investigated by many groups, 
thus far taking one of two approaches, either utilizing a 
transmission line composed of a series array of Joseph¬ 
son junctionsi^AiiiE or a transmission line utilizing the 
nonlinear kinetic inductance of a narrow superconduct¬ 
ing wireiid^i^S. These investigations have revealed the 
need for engineered dispersion to be introduced into the 
transmission line to facilitate phase matchingiiddiidiSi. 
Designs which utilize periodic loadingiid^idd and the ad¬ 
dition of resonant element a^^d^ to facilitate phase match¬ 
ing have shown promise, however at the expense of in¬ 
creased complexity, higher tolerances, and longer propa¬ 
gation lengths (2 cm — 1 m). We propose an alternative 


approach and utilize the nonlinear properties of a one 
dimensional chain of coupled asymmetric SQUIDs as a 
transmission line in a TWPA to achieve phase matching 
and show that exponential gain can be realized over a 
wide bandwidth. 

The proposed TWPA utilizes the tunable nonlinear¬ 
ity of a one dimensional chain of asymmetric SQUIDs 
with nearest neighbor coupling through mutually shared 
Josephson junctions as a transmission line to overcome 
phase matching limitations. A magnetic flux $ threads 
each SQUID to allow for tunability of the linear and non¬ 
linear properties of the transmission line. A weak signal 
to be amplified and a strong pump tone will be incident 
on the input port to the transmission line. Due to the 
nonlinearity of the transmission line the weak signal will 
undergo parametric amplification through a degenerate 
four-wave mixing (FWM) proces o^^i^^ . The amplifica¬ 
tion process is the most efficient when the total phase 
mismatch is close to zero. However due to the nonlin¬ 
earity of the transmission line a strong pump modifies 
phase matching through self and cross phase modulation 
resulting in a phase mismatch. The linear dispersion of 
the transmission line along with spectral separation of 
the signal and pump angular frequencies can be used 
to compensate for the nonlinear phase mismatch. The 
unique feature of the proposed TWPA is that the linear 
and nonlinear dispersion can be tuned with $, and the 
nonlinearity can even change sign. By adjusting $ for a 
given pump power, phase matching can be achieved. 


II. TWPA DESIGN 

The design of the proposed TWPA is shown in Figure 
[T](a). Each cell of the transmission line is an asymmet¬ 
ric SQUID with a single ’’small” Josephson junction with 
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critical current Ij^o and capacitance Cjs in one arm and 
two ’’large” Josephson junctions with critical current IjiQ 
and capacitance Cji in the other arm. 


small junctions. An expansion of the current phase rela¬ 
tion of the current flowing through the backbone of the 
unit cell around (/j = 0 is given by 
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FIG. 1. (color online). TWPA based on a chain of asymmet¬ 
rically coupled SQUIDs (a) Circuit schematic of the proposed 
TWPA. Each unit cell of the TWPA is threaded with a mag¬ 
netic flux $ and has a parasitic capacitance to ground Cgnd- 
The coupled SQUID transmission line can realistically be im¬ 
plemented in a coplanar transmission line geometry. The ge¬ 
ometrical inductance of the transmission line can be made 
negligible in comparison to the Josephson inductance of the 
coupled asymmetric SQUIDs making up the line, (b) Unit 
cell of the TWPA. (c) and (d) Current phase relation of a 
unit cell for various ratios r and 'l>/$o. 


Adjacent cells to one another are coupled via the large 
Josephson junctions. A feature of this arrangement is 
that for an even number of asymmetrical SQUIDs in 
the chain, the Josephson energy remains an even 

function of the phase difference (j) across the chain for 
any value of The proposed TWPA requires a long 
transmission line with many unit cells, this allows to ne¬ 
glect boundary effects and focus on translationally in¬ 
variant solutions. At arbitrary $ the transmission line 
remains symmetric under the translation by two cells. 
The defined unit cell of the transmission line is com¬ 
posed of two large and two small junctions (Figure[IKb)). 
Each unit cell is of length a and has a capacitance to 
ground of Cgnd- The ’’backbone” of the unit cell is 
made up of ’’large” Josephson junctions highlighted in 
red, which are designed to have Josephson energies, EjiQ, 
two orders of magnitude larger than the charging energy, 
Ed = e'^f{2Cji), of the junction. In this case quan¬ 
tum fluctuations of the phase across individual junctions 
are small oc exp[—and thus a classical 
description of the system can be used. The phases on 
the two large junctions for each unit cell are a and a , 
and the total phase across the unit cell \s ip = a + a . 
The Josephson energy Ej{p) of the unit cell is minimized 
when a = a , and a gauge was chosen such that an ex¬ 
ternal magnetic flux would induce phases 27r$/$o on the 



where $o = /i/2e is the flux quantum, h is Planck’s con¬ 
stant, e is the electron charge, and r = Ijio/IjsQ. For 
<i)/$o = 0.5 there exists a critical tq = 4 where the linear 
term in I{p) goes to zero and the cubic term dominates. 
At this point the inductance L oc {dl{p)/dp) of the 
unit cell is the largest. This is the unique feature of the 
superinductor— which is similar by design to the pro¬ 
posed TWPA. As for the proposed TWPA a large induc¬ 
tance of the unit cell is not desired, however tunability 
of the nonlinearity of I{p) is. Fig. [TJc) and[T](d) show 
the current phase relation of a unit cell of the TWPA for 
various r and $. Fig. [TJc) and equation ([IJ show that 
the nonlinearity at I {ip rs 0) is always positive at full 
frustration = 0.5 for r < 16. FigHKd) and equa¬ 

tion © show that for certain r values (example r = 6) 
by adjusting $ the nonlinearity can be tuned over a wide 
range, and can even change sign from negative to posi¬ 
tive. By tuning the nonlinearity it is possible to optimize 
the parametric amplification efficiency of the FWM pro¬ 
cess. 

Using equation o a nonlinear wave equation was de¬ 
rived to describe the node flux pn shown in Fig. ©a) 
along the length of the transmission line^^. Assuming 
constant Cjs, Cji, Cgnd, Ijio, and Ijso along the length 
of the transmission line, no dissipation, and the contin¬ 
uum approximation for a wave-type excitation (A a) 
the following wave equation is derived 


-[- + 2c« 
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( 2 ) 


where L = p^/Ijs^ and p^ = $o/(27r). The first three 
terms of the wave equation represents the linear contribu¬ 
tions to the dispersion on the transmission line due to the 
distributed inductances and capacitances and how they 
can be tuned with r and $. The fourth term describes 
the nonlinearity and how the nonlinear coupling coeffi¬ 
cient, 7 = [a^/((y9oU)] [(r/48) J- (1/3) cos (27r<I>/$o)] can 
be tuned with $. The solution to equation ([2|) is assumed 
to be four traveling waves, where in the degenerate case 
the two pump angular frequencies ujp are equal, a sig¬ 
nal LOs, and a generated idler tone Wi = 2uip — ujg. By 
taking the slowly varying envelope and undepleted pump 
approximations a set of coupled mode equations is de¬ 
rived to describe the propagation of the signal and idler 
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traveling waves: 


das 

dz 

daj 

dz 


ii-^klkiks{2kp - h)a*\Apo\'^ 
Suj'^Cgnd 

i3^k'^ksh{2kp - ks)al\Apo\'^ 
SUJ^Cgnd 


= 0 , 

(3) 

= 0 , 

(4) 


where as and are the complex signal and idler ampli¬ 
tudes, km is the wave vectors of the pump, signal, and 
idler (m = {p, s,i}) (see Appendix A), Afc = ks + ki — 2kp 
is the phase mismatch due to the linear dispersion in the 
transmission line, the total phase mismatch including self 
and cross phase modulation is At = —Ak + 2 q!j, — as — ai, 
and |Apo| is the undepleted pump amplitude. The terms 
am oc 7 |ApoP, are the self and cross phase modulation of 
the wave vectors per unit length. 

Equations m and o are similar to well established 
fiber parametric amplifier theory^i^ and have the fol¬ 
lowing solution to describe the power gain of the signal 
in the presence of a strong pump with zero initial idler 
amplitude Gs = \ cosh( 5 z) —i{Kj2g) sinh(pz)p where the 
exponential gain factor is 


k‘^kf{2kp - ks){2kp - ki)ui^ 
k^uj'iuj^ 

At nearly phase matching conditions k ~ 0, p is posi¬ 
tive and real and the signal gain has an exponential de¬ 
pendence on the length of the transmission line Gs = 
|gg 2 / 2 | 2 _ Pqj. pump amplitudes |Apo| — lA^ol and 

LOs = ujp phase matching is not a concern. For larger 
pump amplitudes the parametric amplification process 
losses phase matching through self-phase modulation of 
the pump, when this happens g becomes small and k 
is large, the gain scales quadratically as a function of 
TWPA length. 

The phase mismatch due to the large pump amplitude 
can be compensated with the linear dispersion along the 
transmission line if ap and Ak are of opposite sign. 
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III. NUMERICAL SIMULATIONS 

From this point on numerical results are presented for a 
realizable set of parameters for the proposed TWPA. For 
each unit cell r = 6 , Ggnd = 50 fF, Gjs = 50 fF, Gji = 
rCjs, Ijso = 1 AiA, and Ijio = rijso- In choosing Ggnd and 
the inductance of the large junctions Lji = $o/(?’Ijso) 
which ultimately sets the characteristic impedance of the 
transmission line special attention was made to achieve 
an impedance near 50 O over the tunable range of the 
TWPA in order to maintain compatibility with commer¬ 
cial electronics. A realizable unit cell size based on our 
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FIG. 2. (color online) Pump tone wave vector and pump self¬ 
phase modulation per unit length a as a function of magnetic 
flux ^/'Fq. The inset shows the change in sign of ap versus 


fabrication process is o = 8 The traveling waves 

used in the numerical results are a pump tone with an¬ 
gular frequency uip/{2 tt) = 6.5 GHz and power —76 dBm 
which is equivalent to Iprms ~ 0.6/iA , the signal angu¬ 
lar frequency Wg was varied in most cases and the idler 
angular frequency is uJi = 2ujp — ujs with initial signal 
and idler power levels 80 dB and 160 dB lower than the 
pump power respectively. 

Shown in Figure [2] is the dependence of kp and ap on 
<i)/<i)o. From the inset in Figure [5] it can be seen that 7 
and as a result ap changes sign from positive to negative 
for d) > 0.31<i)o and more importantly is of opposite sign 
to Afc > 0 (see Appendix A) for this transmission line. 
By adjusting <1) and at a particular cOp and uis it is possible 
to utilize Afc which increases with $ in the transmission 
line to compensate the phase mismatch due to self-phase 
modulation of the pump. 

Fig [31(a) shows numerical simulations of the signal 
gain as a function of signal frequency for the proposed 
TWPA with a transmission line length of 600a. For a 
magnetic flux tuning of $/$o = 0.45 and pump power 
—76 dBm, Fig. (HKa) (red line) there are two regions 
uJs/{2Tr) = 3.6 and 9.4 GHz where perfect phase match¬ 
ing K = 0 can be achieved, and for comparison the phase 
mismatch dependence on signal frequency is shown in 
Fig. IHKb) (red line). For k fv 0 and g is real and 
large, the gain depends exponentially on the TWPA 
length Fig. jSKc) (solid red line). When the phase mis¬ 
match is the largest at uJs/{2'k) = 6.5 GHz, g is small 
and K is large, the gain depends quadratically on the 
length of the TWPA as shown in Fig. |3Kc) (dashed red 
line). Under the phase matching conditions there exist 
regions of exponential gain and quadratic gain depending 
on signal frequency, the 3 dB bandwidth of the TWPA 
^1.5 GHz is limited to two regions where k 0 centered 
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FIG. 3. (color online) Calculated gain of the proposed TWPA. 
In all three panels the color red and blue represent flux tun¬ 
ings (pump powers) of 'I>/4>o = 0.45 (—76 dBm) and <l'/<I>o = 
0.5 (—73 dBm) respectively, (a) Signal gain in dB as a func¬ 
tion of signal frequency, (b) Phase mismatch as a function of 
signal frequency, (c) Dependence of the signal gain as a func¬ 
tion of transmission line length in units a. Solid red and blue 
lines correspond to cUs/(27r) = 9.1 GHz and 8.4 GHz respec¬ 
tively. Dashed red and blue lines correspond to ijJs/{2'k) = 
6.5 GHz. As can be seen the region of exponential gain is 
K ~ 0. According to equation © the most optimal gain does 
not necessarily occur at perfect phase matching k = 0 due 
to the pre-factor to a^. At signihcant k the gain depends 
quadratically on the length shown in panel (c) dashed lines 
and solid blue line. For the flux tuning of 4>/4>o = 0.5 expo¬ 
nential gain is impossible at all frequencies. 


at uJs/{2tt) = 3.9 and 9.1 GHz.With magnetic flux tun¬ 
ing value of <f>/$o = 0.5 and pump power —70 dBm, k 
is large Fig. |3Kb) (blue line) and only a quadratic gain 
dependence is possible. Since the signal gain increases 
quadratically for all frequencies a relatively flat gain char¬ 
acteristic of the amplifier can be achieved with a sig¬ 
nal gain of 23 dB over a 3 dB bandwidth greater than 
5.4 GHz. 

The saturation power of the TWPA is limited by pump 
depletion effects, which generally occurs when the sig¬ 
nal amplitude becomes comparable to the pump. At 
such amplitudes the pump is depleted and the gain of 
the TWPA decreases. In order to evaluate these effects 



FIG. 4. (color online) Calculated 1 dB compression point 
and minimum TWPA length as a function of magnetic flux 
<l?/4>o to maintain a signal gain of 20 dB. As ‘F/'Fo varies 
from 0.35 to 0.5 the pump power was varied from —68 dBm 
to —76 dBm to maintain phase matching conditions between 
a signal and pump tone at angular frequencies uJs/{2n) = 
6 GHz and uip/{2'K) = 9 GHz respectively. The inset shows 
the signal gain as a function of signal power for 4>/<l>o = 0.37 
and device length 341a, the 1 dB compression point occurs at 
-98 dBm. 


the coupled mode equations (see Appendix B) without 
the un-depleted pump approximation and taking into ac¬ 
count self- and cross-phase modulation are solved to de¬ 
termine the real amplitude and phase mismatch as a func¬ 
tion of z along the length of the transmission line. Figure 
m (inset) shows how the signal gain decreases with the sig¬ 
nal power due to pump depletion effects, for = 0.37 
and a TWPA length of 341o. The signal gain and 
the phase mismatch depend on magnetic flux through 
Qfp oc qfcpIApoP and km. For each <i)/<i)o the pump power 
is varied over the range —76 dBm {Iprms ~ 0.1/j/o) to 
—68 dBm {Iprms ~ 0.25/j7o) to maintain phase matching 
at ojs/ (27r) = 6 GHz with a uip/{2Tr) = 9 GHz pump. For 
each <i)/<i)o the minimum length of the TWPA to achieve 
a signal gain of 20 dB is shown in Fig. |4](red line). When 
ap and km decrease with decreasing $/$o a stronger 
pump is required to maintain phase matching conditions 
which results in a larger PidB up to —98 dBm, smaller 7 , 
and a longer transmission line to maintain a signal gain of 
20 dB. When ap and km increases with <i)/<i)o a smaller 
pump power is required to achieve phase matching, PidB 
decreases, 7 increases resulting in a shorter transmission 
line with a minimum length of 211 a to achieve a signal 
gain of 20 dB. 


IV. SUMMARY 


In conclusion a TWPA design based on a chain of cou¬ 
pled asymmetric SQUIDs has been presented. The pro¬ 
posed design allows for great flexibility where a magnetic 
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flux can be used to tune the nonlinearity of the trans¬ 
mission line to achieve phase matching conditions in a 
four-wave mixing process. Numerical simulations have 
shown that the proposed amplifier can achieve gains of 
23 dB with a 3 dB bandwidth of greater than 5.4 GHz 
at a center frequency of 6.5 GHz. Under different tuning 
conditions gains of greater than 20 dB can be achieved 
with a minimum transmission line length of 211 unit cells 
and a saturation power of up to —98 dBm with 341 unit 
cells. The proposed amplifier is ideally suited for mul¬ 
tiplexed readout of quantum bits or kinetic inductance 
based astronomical detectors. 
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Appendix A: Analytical Approximation 


In this section we derive an analytical solution to the 
coupled nonlinear wave equations used to describe the 
interaction between the signal and pump traveling waves 
propagating along the transmission line composed of a 
chain of coupled asymmetric SQUIDs. The energy phase 
relation of a unit cell Figure 1(a) is 


Ej{ip) = -2EjiQCOS 


Ejs^s cos 


(^) -Ej,ocos 

(---a. 



27r<i)\ 


(Al) 


Expanding dEj{(p)/dip\^=o gives the approximate cur¬ 
rent phase relation describing the current Inif) flowing 
through the backbone of unit cell n Eq. (1). Utilizing 
Kirchhoff’s current law In-i{(p) — Cgnd dP‘f / dt^ = Inif) 
and neglecting the effects of dissipation and assuming 
a sufficiently long wavelength A ^ a of the signal and 
pump the following wave equation is derived for a posi¬ 
tion z along the transmission line 


2 

rr 

Ti2 


+ 2 cos 





a Cjs 



dPdz^ 




C'g„d ^^2 ^'dzi\dzJ 




= 0 , 


(A2) 


where L = ifo/Ijs, Vo = ‘&o/27r, r = Ijio/Ijso, and 
7 = [a'^/(i/?§L)]|(r/48) -|- (1/3) cos (27r<I>/$o)]. The so¬ 
lution to Ea. (lA2l) is assumed to be a superposition of 
a pump, signal, and idler traveling waves propagating 


along the transmission line of the form 


+C.C 


(A3) 


where c.c. denotes complex conjugate, Am is the complex 
amplitudes, km is the wave vectors, and uim is the angular 
frequencies of the pump, signal, and idler (m = {p, s, i}). 
A degenerate four-wave mixing process is considered un¬ 
der the following frequency matching condition uis+coi = 
2ujp. Ea. (IA3p is substituted into Ea. (jA2ll and assuming a 
slowly varying envelope of the propagating waves where 
\d'^Am/dz'^\ <C \kmdAm/dz\ and \dAm/dz\ < \kmAm\ 
and a uniform transmission line where Cgnd, Cjs and km 
are constant, a set of coupled mode equations which de¬ 
scribes the propagation of the pump, signal, and idler 
waves along the transmission line is determined: 


(A4) 


r) A 

- lapAp = 0 , 

dAs ^ '^‘d‘'ykpkiks{2kp — ki)A^Ap ^ 

XCXs-^s C — U, 


dz 

dA, 

dz 


xoijAi 


Scd^sCgnd 


i2,-fklkski{2kp - ks)A*Al 


(A5) 


^U^iCgnd 


= 0 , 


(A6) 


where a large pump amplitude relative to the signal and 
idler amplitudes was assumed, decoupling the pump, and 
the quadratic terms in Agj were neglected, Afc = fcs -b 
ki — 2kp is the phase mismatch due to linear dispersion, 
and am is the self-phase modulation per unit length a: 


'iCgndOJs 


^Cgnd^f 


(A7) 


37^pl^poP 

SCgr^dUJ^ ’ 

where |Apo| is the initial pump amplitude. The linear 
dispersion relation for this transmission line is 

km = U;mVLCg„d ^ 

«V[i+2cos(27r^)]-<LQ«(i + 2) 

Assuming an un-depleted pump amplitude and the fol¬ 
lowing substitutions Ap{z) = Apoe*“r^ solution to Eq. 
(IA4I) . As{z) = as{z)e^°‘^^, and Ai{z) = ai{z)e^°'^ into 
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(IA5I) and (IA6I) to obtain: 

das ii'lk'^kiks{2kp - h)a*\Apo\'^ _ 

dz BujlCgnd ^ 


(A9) 


dai ^‘^^k^kski{2kp A:s)Q,g|Apo| 
dz SujfCgnd ^ 


(AlO) 


where k = —Ak + 2ap — is the total phase mis¬ 

match. Equations (IA9I) and (lAlOp are similar to well 
established fiber parametric amplifier theory and have 
the following solution to describe the amplitude of the 
signal along the length of the transmission line assuming 
zero initial idler amplitude: 


as{z) 




cosh(5z) 


— sinh(gz) 


giKz/2 


(All) 


A similar solution to (lAlip exists for the idler amplitude. 
The exponential gain factor is 


g = 


I fk‘jk^{2kp - ks){2kp - ki)uj^ 


k^uj^ujs 


- '“p - 


(A12) 

The signal power gain can be determined from Eq. lAllI 


Appendix B: Numerical Analysis 


Below the solutions to the coupled mode equations 
which govern the degenerate four wave mixing process 
between the signal, idler, and pump tones and pump de¬ 
pletion effects are presented. Plugging Eq. IA3I into Eq. 
1X2] and making the slowly varying envelope approxima¬ 
tion the following coupled mode equations are derived: 
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cosh( 5 z)-sinh(gz) 

2ff 


(A13) 


Eor the proposed chain of asymmetric SQUIDs the phase 
mismatch due to linear dispersion is always real and non¬ 
negative for r > tq. Using the frequency matching con¬ 
dition and Eq. IA8I we show that 


Ak = Cl 


CUs 




-f 


LOi 


C2Wi 




C2W" 


where ci 


V ^^qnd , 

- and 

ayi + 2cos(27r^) 



LCgs{^ + 2) 

(i + 2cos(27r|^))’ 


Ak Ri 3ciC2Aw^Wp > 0, 


where Aw = uis—ujp = Wp—uji from the frequency match¬ 
ing condition. We also show that the sign of the first term 
in Eq.(6) is only dependent on ap for r > rg, 


kL 






(1 - C2UJm) '' 






w2 ^ w2 uj 


ClWpf 1 -I- 9c2Aw^ 



(A14) 


> 0 . 


The coupled complex differential equations are solved 
by converting the complex amplitudes to Am = 
Bm{z)e^^"‘G) and finding the solutions for the real ampli¬ 
tudes Bp{z), Bs{z), Bi{z), and the total phase mismatch 
0s(z)+9i(z)—29p(z)+Ak using the Runge-Kutta method. 


Appendix C: Parameter Variations 

In this section we determine how variations in Joseph- 
son junction parameters and tuning magnetic flux gradi¬ 
ents effect the operation of the proposed TWPA. Both 
types of variations manifest themselves as a position de¬ 
pendent alteration to both the linear and nonlinear dis¬ 
persion from unit cell to unit cell which ultimately effects 
the phase matching and the signal gain in the TWPA. 
Since the transmission line of the TWPA is made up 
of coupled asymmetric SQUIDs variations in one unit 
cell through mutual coupling, effects the superconduct¬ 
ing phase distribution in neighboring unit cells. 

We model the variations in Josephson junction param¬ 
eters as a normal distribution with mean critical currents 
Ijso and Ijio for small and large Josephson junctions re¬ 
spectively, and a relative standard deviation ajj which 
applies to both small and large junctions. The proposed 
TWPA was modeled by minimizing the total Josephson 
energy of a 200 unit cell chain by adjusting the phases on 
all the large junctions. A gauge was chosen such that the 
external magnetic flux would induce phases 27r<I>/<I>o on 
the small junctions. The constraint of the minimization 
procedure was that the sum of the phase drops on all 
of the large junctions would be equal to the total phase 
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across the chain. From the minimization a position de¬ 
pendent Ak„i and A 7 from unit cell to unit cell was de¬ 
termined, vectorized, and introduced into the numerical 
analysis of the coupled mode equations. 




FIG. Al. Simulated signal gain as a function of the relative 
standard deviation of the critical currents of both the small 
and large Josephson junctions which make up the TWPA. 
Simulations were performed at "F/iho = 0.45 (blue circles) 
and 4>/<I>o = 0.5 (red squares), u}p/{2n) = 9 GHz, a;s/(27r) = 
6 GHz, TWPA length 200a, and r = 5 and r — 6 for Panel (a) 
and (b) respectively. The error bars represent the variation 
in signal gain due to the stochastic solution of the coupled 
mode equations to a normal distribution of junction parame¬ 
ters along the length of the TWPA. 

Numerical simulations of the TWPA which take into 
account a normal distribution in Josephson junction pa¬ 
rameters were performed on the shortest 200a TWPA 
with a uJp/{2tt) = 9 GHz and pump power of —76 dBm. 
Shown in FIG. lAH a) and FIG. lAlI b) is the signal gain 
at uJs/{2Tr) = 6 GHz) where k w 0 as a function of 
the relative standard deviation ajj for different r and 
<i>. For each data point 50 numerical simulations were 
performed and the error bars represent the spread in 
signal gain due to the stochastic nature of the solu¬ 
tion to the coupled mode equations to a normal distri¬ 
bution in junction parameters. FIG. lAlI a) shows for 
r = 5 and = 0-5 the tolerance of the proposed 

TWPA to a variation in junction parameters is limited 



mOg per TWPA unit cell 


FIG. A2. Simulated signal gain as a function of magnetic 
flux gradient per unit cell across a TWPA composed of 200 
unit cells, for magnetic flux tunings of cO/iho = 0.45 (blue 
circles) and "F/iho = 0.5 (red squares). Error bars represent 
the stochastic solutions of the coupled mode equations due to 
a normal distribution in Josephson junction parameters with 
ajj = 2.4%. 



Signal Frequency (GHz) 


FIG. A3. Numerical simulation of the signal gain as a 
function of signal frequency for pump frequencies uip/ (27r) = 
3 GHz (green line), 5 GHz (red line), and 9 GHz (blue line) 
respectively, pump power of —76 dBm, r = 6, TWPA length 
of 600a and 'I>/4>o = 0.45. 


to ajj = 1%, where the signal gain drops by 1 dB. FIG. 
lAlI b) shows for r — 6 and <I>/<I>o = 0-5, the tolerance in¬ 
creases to ajj = 3% which can be realized with present 
day fabrication technology where the on-chip variation 
in junction parameters follows a normal distribution and 
tj = 2.4 — . When the magnetic flux is tuned 




















to $/$o = 0.45 the nonlinearity is decreased, however, 
the tolerated ajj increases to 4.5%. A trade-off exists 
between the magnitude of the nonlinearity (for small r 
or large $) and the tolerance of the TWPA to a variation 
in junction parameters. 

Numerical simulations of the signal gain as a func¬ 
tion of magnetic flux gradient per unit cell across the 
TWPA is shown in FIG. IA2l for tuning magnetic fluxes of 
<i)/<i)o = 0.45 (blue circles) and = 0.5 (red squares). 
The error bars represent the stochastic nature of the nu¬ 
merical simulations of the coupled mode equations due 
to a normal distribution in Josephson junction parame¬ 
ters with ajj = 2.4%. At 4>/$o = 0.5 the TWPA has 
a 1 dB drop in signal gain for a magnetic flux gradient 
of 1 m<i>o per unit cell. When the magnetic flux is tuned 
to = 0.45 where the TWPA is more sensitive to 

magnetic flux variations, the signal gain drops by 1 dB 
at 0.5 m$o per unit cell. In our measurement setup with 
moderate magnetic shielding we have observed magnetic 


flux gradients of 24 ^Gauss/ — which is equiva¬ 
lent to 0.1 m<i)o per unit cell of the TWPA. 

Appendix D: Different Pumping Conditions 

Numerical simulations of the signal gain as a function 
of signal frequency is shown in FIG. IA3I for alternative 
pump frequencies Up/ (27r) = 3 GHz, 5 GHz, and 9 GHz 
to be compared with the results shown in FIG. 3(a). Un¬ 
der phase matching conditions, k Rs 0, the maximum in 
the signal gain depends exponentially on g (see Eq. 5), 
which has a strong dependence on ujp. FIG. IA3I shows 
that the maximum gain and bandwidth increase with ujp 
and thus shorter transmission lines can be realized as was 
demonstrated in FIG. 4. For smaller ujp the maximum 
gain and bandwidth decreases and thus longer transmis¬ 
sion lines are required even under optimal tuning condi¬ 
tions for such pump frequencies. 
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